Qubit-oscillator dynamics in the dispersive regime: 
analytical theory beyond rotating-wave approximation 
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We generalize the dispersive theory of the Jaynes-Cummings model beyond the frequently em- 
ployed rotating-wave approximation (RWA) in the coupling between the two-level system and the 
resonator. For a detuning sufficiently larger than the qubit-oscillator coupling, we diagonalize the 
non-RWA Hamiltonian and discuss the differences to the known RWA results. Our results extend 
the regime in which dispersive qubit readout is possible. If several qubits are coupled to one res- 
onator, an effective qubit-qubit interaction of Ising type emerges, whereas RWA leads to isotropic 
interaction. This impacts on the entanglement characteristics of the qubits. 

PACS numbers: 03.67.Lx, 03.70.+k, 42.50.Hz, 42.50.Pq 
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I. INTRODUCTION 

More than forty years ago, Jaynes and Cummings [1| 
introduced a fully quantum mechanical model for the in- 
teraction of light and matter, which are represented by a 
single harmonic oscillator and a two-level system, respec- 
tively. Within dipole approximation for the interaction, 
that model is expressed by the Hamiltonian 

H =^-a z + hLoa)a + hga x (a) + a), (1) 

where he is the level splitting of the two-level system, 
henceforth "qubit" , uj is the frequency of the electromag- 
netic field mode, and g the dipole interaction strength. 
The Pauli matrices cr Q , a = x,y,z, refer to the two- 
level system, while <v and a denote the bosonic creation 
and annihilation operators of the electromagnetic field 
mode. This model describes a wealth of physical phe- 
nomena rather well and by now is a "standard model" 
of quantum optics. A particular experimental realiza- 
tion of the Hamiltonian is an atom interacting with the 
field inside an optical cavity, usually referred to as cavity 
quantum electrodynamics (cavity QED). Corresponding 
experiments have demonstrated quantum coherence be- 
tween light and matter manifest in phenomena such as 
Rabi oscillations and entanglement [2], |3| ■ 

Despite its simplicity, the Hamiltonian |T]) cannot be 
diagonalized exactly and, thus, is often simplified by 
a rotating-wave approximation (RWA). There, one ex- 
presses the qubit-cavity interaction in terms of the ladder 
operators cr = \{o x ± ia y ). In the interaction picture 
with respect to the uncoupled Hamiltonian, the coupling 
operators a + a, a~a^ and a~a, a + a) oscillate with the 
phase factors exp[±i(oj — e)t] and exp[±i(u + e)t] , respec- 
tively. Operating at or near resonance, the cavity-qubit 
detuning is small, |e — w| <C e + so that the former op- 
erators oscillate slowly, whereas the latter exhibit fast 
"counter-rotating" oscillations. If in addition, the cou- 
pling is sufficiently weak, g <min{£,w}, one can sepa- 
rate time scales and replace the counter-rotating terms 



by their vanishing time average. Then one obtains the 
Jaynes-Cummings Hamiltonian [l[ 

Hjxwa — ~~2° 7 ' + hiua'a + hg(a~ a) + a + a) . (2) 



Lately, new interest in Jaynes-Cumming physics has 
emerged in the solid state realm. There, one implements 
artificial atoms with Cooper-pair boxes (charge qubits) 
[H or superconducting loops (flux qubits) [$\. The role 
of the cavity is played now by a transmission line or 
a SQUID depending on the architecture @, 0] , or even 
a nanomechanical oscillator J8|. Since the first experi- 
mental realizations in 2004 0, @ , a plethora of results 
has been obtained, such as quantum-non-demolition-likc 
readout of a qubit state Q , the generation of Fock states 
[ir]| |. the observation of Berry phases multi-photon 
resonances M , entanglement between two qubits inside 
one cavity 0^ , [l3[ , and the demonstration of a two-qubit 
algorithm [lj . 

These experiments have in common that they operate 
in the strong coupling limit, that is, the coupling g is 
larger than the linewidth of the resonator. On the other 
hand, g is typically two orders of magnitude less than 
the qubit and resonator frequencies. In this scenario the 
Jaynes-Cummings model @ has been shown to describe 
the experiments faithfully. 

Of practical interest is the dispersive limit, in which 
the qubit and the resonator are far detuned compared to 
the coupling strength, g<c|e — u)\ [1, [IB]. In this regime 
a non-demolition type measurement of the qubit can be 
performed by probing the resonator [J, fl6| . Moreover, 
it is possible to simulate quantum spin chains with two 
or more qubits that are dispersively coupled to one res- 
onator [lj, [3, G3]- The complementary architecture 
of two cavities dispersively coupled to one qubit allows 
building a quantum switch [18|. All these ideas have 
been developed from the RWA model ([2]) in the dispersive 
limit or from according generalizations to many qubits or 
many oscillators. Thus, these theories are restricted to 
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the range 

.9 < \e-w\ <e + w, (3) 

where the first inequality refers to the dispersive limit, 
while the second one has been used to derive the RWA 
Hamiltonian {2} from the original model 

In recent experiments, efforts are made to reach an 
even stronger qubit-cavity coupling g. Thus, it will even- 
tually be no longer possible to fulfill both inequalities 
([3]) [19 ]. In particular, when trying to operate in the 
dispersive limit, the second inequality may be violated, 
so that RWA is no longer applicable. Non-RWA effects 
of the model Hamiltonian (JTJ) have already been studied 
in the complementary adiabatic limits e <^uj and uj -C e 
[20I . I2H l22l | . Furthermore, Van Vleck perturbation theory 
has been used in the resonant and close-to-resonant cases 
[23|]. Finally, polaron transformation (24|, cluster meth- 
ods [iH , wave-packet approach [2H ] and even generalized 
RWA approximations |27| have been considered. 

Motivated by the importance of the dispersive regime 
and in view of the experimental tendency towards 
stronger qubit-oscillator coupling, we present in this work 
a dispersive theory beyond RWA, so that the second con- 
dition in Eq. ([3]) can be dropped. This means that our 
approach is valid under the less stringent condition 

ff«|e-w|, (4) 

which implies that the detuning is not necessarily smaller 
than e and u>. In order to set the stage, we briefly review 
in Sec. [TT] the dispersive theory within RWA. In Sec. IIII1 
we derive a dispersive theory for Hamiltonian ([T]) beyond 
RWA, which we generalize in Sec. IIVI to the presence of 
several qubits. 



A [-ffRWA j A_ ] + 1 A 2 [ [i?RWA , A_ ] , A_ ] , which can be eval- 
uated to read 

He , hq 2 „ ( hq 2 S\ + 
#RWA,disp = ycr + ^-cr z + I Huj + -j^-<r z J a T a . (9) 

The physical interpretation of is that the oscillator 
frequency is shifted as 

uj^uj±g 2 /A, (10) 

where the sign depends on the state of the qubit. If 
one now probes the resonator with a microwave signal 
at its bare resonance frequency u>, the phase of the re- 
flected signal possesses a shift that depends on the qubit 
state. This allows one to measure the low-frequency dy- 
namics of the qubit @. Since, according to Eq. (|9|), the 
qubit Hamiltonian (he/2)a z commutes with the disper- 
sive coupling (Hg 2 / A)a z a f a, this constitutes a quantum 
non-demolition measurement of the qubit,which has al- 
ready been implemented experimentally 0, H2, HH, H^j • 
In turn, the qubit energy splitting is shifted depending 
on the mean photon number n — (a^a). Accordingly, one 
can also measure the mean photon number, and even per- 
form a full quantum state tomography of the oscillator 
state [13]. Note also that besides the condition of A be- 
ing small, the perturbational result ([9]) is accurate only if 
the mean photon number (n) does not exceed the critical 
value ricrit = 1/4A 2 . For larger photon numbers, higher 
powers of the number operator a^a must be taken into ac- 
count [3~il [HJ • Henceforth, we restrict ourselves to the 
so-called linear dispersive regime, in which the photon 
number is clearly below the critical value ri cr i t . 

III. DISPERSIVE REGIME BEYOND RWA 



II. DISPERSIVE THEORY WITHIN RWA 

The dispersive limit is characterized by a large detun- 
ing A = e— uj as compared to the qubit-oscillator coupling 
g. Thus, 



(5) 



represents a small parameter, while the RWA Hamilto- 
nian @ is valid for |e — u)\ <C e + uj. Then it is convenient 
to separate the coupling term from the RWA Hamilto- 



nian, i.e., to write -Hrwa 
butions 

Hq = ~^° Z + hjja) a , 
X± =<j~a) ± a + a . 
Applying the unitary transformation, 



Hq + hgX + with the contri- 
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RWA 



(0) 
(7) 

(8) 



one obtains for the transformed Hamiltonian -ffdisp 



r)t H 



RWA 1 1 RWA ^ RWA 



to second order in A: H, 



disp 



H 



RWA- 



It is now our aim to treat the original Hamiltonian ([T]) 
in the dispersive limit accordingly, i.e., to derive an ex- 
pression that corresponds to Eq. §§§ but is valid in the 
full dispersive regime defined by inequality Q. Going 
beyond RWA, we have to keep the counter-rotating cou- 
pling terms 



Y+ = a+a 1 ± cr" 



(11) 



which are relevant if either of the relations g <C min{e, lo} 
or e — lo\ <C e + lo is violated. Separating again the 
qubit-oscillator coupling from the bare terms, we rewrite 
Hamiltonian |T]) as 



H = H + hgX+ + HgY + 



(12) 



which differs from -Wrwa by the last term. It will turn out 
that a unitary transformation corresponding to Eq. (J5]) 
is achieved by the operator 
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_ „AX_+AY_ 



Here we have introduced the parameter 
9 9 



A 



2e - A ' 



(13) 



(14) 
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which obviously fulfills the relation A < A, since e and 
co are positive. Thus, whenever A is small, A is small as 
well. Nevertheless, under condition ((4|), A and A may be 
of the same order. 

Proceeding as in Sec. HH we define the dispersive 
Hamiltonian i?disp = HD. Using the commutation re- 
lations [Y+,Y_] = o- z (2at a + l)-l, [e/2<7 z + wa + a,Y_] = 
-(e + u)Y + and \Y±,X T ] = a z {a 2 + (at) 2 }, we obtain 
the expression 



fr e z t- t 

-Hdisp = ycr + hioa [ a 
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(15) 



which is valid up to second order in the dimensionless 
coupling parameters A and A. 

As compared to the RWA result (jSJ), we find two dif- 
ferences: First, the prefactor of the coupling has a con- 
tribution that obviously stems from A. Second and more 
importantly, the coupling is no longer proportional to the 
number operator a^a, but rather to (a' + a) 2 . Thus, the 
operator Y± has turned into the counter-rotating con- 
tributions a 2 and (a') 2 . For this reason, the dispersive 
Hamiltonian ([IB"]) is not diagonal in the eigenbasis of the 
uncoupled Hamiltonian Hq. 

Nevertheless, it is possible to interpret the result as 
a qubit-state dependent frequency shift by the following 
reasoning. Let us interpret Hcoa'a as the Hamiltonian of 
a particle with unit mass in the potential \co 2 x 2 , where 
x = ■\/h/2uj(a' + a). Then the qubit-oscillator coupling 
in Eq. (|T5"|) modifies the potential curvature to 2 , such that 
the oscillator frequency undergoes a shift according to 
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FIG. 1: Oscillator frequency shift as function or the qubit 
splitting e = uo + A for the spin state |j) obtained (a) within 
RWA, Eq. dTDJ, and (b) beyond RWA, Eq. (TF]|. The lines 
mark the analytical results, while the symbols refer to the 
numerically obtained splitting between the ground state and 
the first excited state in the subspace of the qubit state |J.}. 



1± 



2ff 2 



CO 



2e- A 



(16) 



Again the sign depends on the qubit state. To be con- 
sistent with the second-order approximation in g, we 
have to expand also the square root to that order. This 
complies with the experimentally interesting parameter 
regime where g < co. We finally obtain 



CO 



io±g 2 



1 



2e- A 



(17) 



As for the RWA Hamiltonian, we find that the qubit state 
shifts the resonance frequency of the oscillator. This re- 
sult is not only of appealing simplicity, but also has a 
rather important consequence: Dispersive readout is pos- 
sible even when the qubit-oscillator coupling is so strong 
that condition cannot be fulfilled, that is, when the 
RWA result is not valid. 

For a quantitative analysis of our analytical findings, 
we compare the frequency shifts (fTO]) and (fT7)) with nu- 
merical results. In doing so, we diagonalize the Hamil- 
tonian |T|) in the subspace of the qubit state ||), where 
<J z \l) = — 14). The results are depicted in Fig. [1] 



For a qubit splitting e close to the cavity frequency to, 
i.e., outside the dispersive regime, the analytically ob- 
tained frequency shifts diverge. This behavior is cer- 
tainly expected for an expansion in <?/A. For a rela- 
tively small coupling g/co< 0.025, the RWA result [panel 
(a)] agrees very well with the numerical data in the dis- 
persive regime. In the case of larger coupling strengths, 
g/co > 0.05, the predictions from RWA exhibit clear differ- 
ences. The general tendency is that RWA overestimates 
the frequency shift for blue detuning A = e — co < 0, while 
it predicts a too small shift for red detuning. 

The data shown in panel (b) demonstrates that the 
treatment beyond RWA yields the correct frequency shift 
in the entire dispersive regime, i.e., whenever the detun- 
ing significantly exceeds the coupling, |A| 3> g. Thus, 
as long as the coupling remains much smaller than the 
oscillator frequency, g -c co, it is always possible to tune 
the qubit splitting e into a regime in which U is fulfilled. 
Moreover, the excellent quantitative agreement of our an- 
alytical result (|17[) with the numerically exact solution in- 
dicates the feasibility to determine g from measurements 
in the strong-coupling limit . 

A particular limit is A — > —co, which corresponds to a 
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vanishing qubit splitting, e — > 0. In this case it is obvi- 
ous from Hamiltonian ((T|) that the coupling to the qubit 
merely entails a linear displacement of the oscillator coor- 
dinate, while the oscillator frequency remains unaffected. 
This limit is perfectly reproduced by our non-RWA result 
(fT7|) . irrespective of the coupling strength. The RWA 
result, by contrast, predicts a spurious frequency shift, 
indicating the failure of RWA. 



IV. SEVERAL QUBITS IN A CAVITY 

An experimentally relevant generalization of the model 
(HJ) is the case of several qubits coupling to the same 
oscillator. The corresponding Hamiltonian reads [33[ 



on a collective coordinate of all qubits. Consequently, 
the cavity response experiences a phase shift from the 
ingoing signal, which in turn contains information about 
that collective qubit coordinate. 



B. Dispersive theory beyond RWA 

As in Sec. lIIII for the one-qubit case, we now extend the 
dispersive theory of the Tavis-Cummings model beyond 
RWA, taking into account the counter-rotating terms of 
the Hamiltonian (I18|) . In analogy to transformation (|13| . 
we employ the ansatz 



D 



>:-\ A • A V 



(22) 



H 



3 



hiua^a + hY^ g^Aa^ + a), (18) 
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where the index j labels the qubits. As for the one- 
qubit case, the rotating wave-approximation is frequently 
applied and yields the Tavis-Cummings Hamiltonian [34j | 

H = h^^<T? + huJa}a + hJ29j x + » ( 19 ) 
where X J ± = (J J a] ± c^o,-, cf. Eq. ©. 



A. Dispersive theory within RWA 



We obtain for each qubit the dimcnsionlcss coupling 
parameter Aj — gj/(cj — to). The dispersive limit is 
now determined by \Xj\ <C 1 for all j. Effective decou- 
pling of the qubits and the cavity to second order is then 
achieved via a transformation with the unitary operator 
exp(— Y^,j XjX 3 _), cf. Eq. ©. The resulting dispersive 
Hamiltonian reads [331 



ffdisp = ^a f a + | ( e i + IT. ) a i 



X! ~r" < " 7 ' 1 Yl J 3k(°jv£ + °-t a k)- 



j>k 



(20) 



Remarkably, the oscillator entails an effective coupling 
between the qubits with the strength 



Jjk = g 3 gk 



l 

A? 



i 



(21) 



which has already been observed experimentally [12j. It 
has been proposed to employ this interaction for building 
qubit networks 1171 and for generating qubit-qubit entan- 
glement [35l l36l. l37l]. Moreover, quantum tomography of 
a two-qubit state has been implemented by probing the 
cavity at its bare resonance frequency [38j |. In this sce- 
nario the oscillator frequency exhibits a shift depending 



where Y 3 _ = aj a — < jjjv and Xj = g/(2e — Aj). Follow- 
ing the lines of Sec. IIII1 i.e. expanding the transformed 
Hamiltonian to second order in A and A, we obtain the 
dispersive Hamiltonian 

H disp = D^HD 

= huja' i a+^Y^ e 3< T * 

3 



3 



+ E J ^K- ( 23 ) 
j>k 

We have introduced the modified coupling strength 



J 



jk 



g,g k \ x - 



1 



1 



2e - A., 



2e - A fc 



> (24) 



which describes the effective interaction between qubits 
j and k, and represents the extension of Eq. (f2"Tj) be- 
yond RWA. The dispersive shifts of the qubit and cav- 
ity frequencies, given by the second and third term of 
Eq. (|2"3"|) , are equally modified as compared to the RWA 
result PD|) . 

Interestingly enough, the effective qubit-qubit interac- 
tion in Eq. ([23]) is of the Ising type cjcr£, whereas RWA 
predicts the isotropic XY interaction cr^CjT + crjcr^ , see 
Eq. (|2"D|) . Thus, the treatment beyond RWA predicts a 
qualitatively different effective model and not merely a 
renormalization of parameters. The Ising term even per- 
sists in the limit 1/Aj 3> l/(2e — Aj). Nevertheless, one 
can recover the RWA Hamiltonian ([20| by writing the 
interaction term as crj 0% — cn + crfcr^ + h.c. and 
performing a RWA for the Ising coupling. This corre- 
sponds to discarding small-weighted, rapidly oscillating 
terms of the type JjkUj'a^ + h.c. 

The difference between the effective models (|2TJ)) and 
([23]) has some physically relevant consequences. First, 
in contrast to the RWA result (|2"0)). Hamiltonian J2"5j) 
does not conserve the number of qubit excitations, which 
will affect the design of two-qubit gates (33[. Moreover, 
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both models possess different spectra, which influences 
entanglement creation. For instance, the ground state of 
the Hamiltonian (|23|) for two degenerate qubits (ei =e 2 ) 
coupled to one cavity is |0)[||J,) — (J/2e)|TT)] an d thus 
exhibits qubit-qubit entanglement. By contrast, the cor- 
responding ground state of the multi-qubit RWA Hamil- 
tonian (f2T))) is the product state |0)||J.). For the case 
of the Hamiltonian (|23p . thermal qubit-qubit entangle- 
ment will consequently be present at zero temperature 
and even at thermal equilibrium [3^ . [40| . 

V. SUMMARY 

We have generalized the dispersive theory for a qubit 
coupled to a harmonic oscillator to the case of far 
detuning. In this limit, it is no longer possible to 
treat the qubit-oscillator interaction Hamiltonian within 
the rotating-wave approximation. Therefore, previous 
derivations need some refinement. It has turned out that 
diagonalizing the Hamiltonian analytically up to second 
order in the coupling constant is possible as well be- 
yond RWA. The central result is that as within RWA, 
the oscillator experiences a shift of its resonance fre- 
quency, the sign of the shift depending on the qubit state. 
In this respect, the difference between both approaches 
seems to be merely quantitative. Nevertheless, our result 
implies an important fact for currently devised qubit- 



oscillator experiments with ultra-strong cavity-qubit cou- 
pling: Dispersive qubit readout is possible as well in that 
regime. The comparison with numerical results has con- 
firmed that our approach is quantitatively satisfactory in 
the whole dispersive regime. 

The corresponding treatment of many qubits coupled 
to the same oscillator is equally possible. In such archi- 
tectures, the oscillator mediates an effective qubit-qubit 
interaction which may be used for gate operations and 
entanglement creation. We have revealed that the form 
of the effective interaction depends on whether or not 
one employs RWA. While RWA predicts an isotropic XY 
interaction, the inclusion of the counter-rotating terms 
yields an interaction of Ising type. This difference im- 
pacts on various proposed entanglement creation proto- 
cols as soon as they operate in the far-detuned dispersive 
regime. 
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